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We briefly review the formulation of chiral quark soliton model and 
explain the difference and similarities with the Skyrme model. Next, we 
apply the model to calculate non-exotic and exotic mass spectra. We con- 
. centrate on large N c counting both for mass splittings and decay widths. 

It is shown that pure large N c arguments do not explain the small width 
of exotic pentaquark states. 
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1. Introduction 
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^ \ There is still a lack of consensus whether the lightest member of the ex- 

' otic antidecuplet has been discovered . Four months after this conference 

results from high statistics Gil experiment at CLAS have been presented 
at the APS meeting with negative result for the photoproduction of <9 + on 
proton 2 . Even more problematic is the sighting of the heaviest members of 
10 that were seen only by NA49 experiment at CERN These states were 
predicted within the chiral soliton models [1J 03 H3 [JJ |Sj . Early estimates of 
the antidecuplet octet splitting, AMjq_ s ~ 600 MeV, obtained in a specific 
modification of the Skyrme model can be found in Ref. jH]. The estimates 
of both + and Stq masses from the second order mass formulae obtained 
in the Skyrme model in 1987 are in a surprising agreement with present ex- 
perimental findings Already at that time, however, the doubts whether 
these predictions were trustworthy had been raised E3 EI • Today they 
were scrutinized and rephrased by other authors |121 1131 114j . 

In 1997 the masses, as well as the decay widths of the exotic states were 
estimated within the chiral quark soliton model 7_. One of the most striking 
predictions of this seminal paper by Diakonov, Petrov and Polyakov [7j was 



(1) 



2 



EPIPHANY PRINTED ON FEBRUARY 2, 2008 



the small width of antidecuplet states. Despite some misprints in this paper 
(see e.g. HE2j) and the model dependent corrections, the narrow width is 
one of the key features of the chiral model predictions which is in line with 
recent experimental findings. 

In this paper we examine the successes and problems of chiral soliton 
models ability to predict properties of exotic baryons. In Section |2] we ar- 
gue that soliton models are in fact quark models and explain the difference 
between quark-soliton models and Skyrme model. Then in Sect. |3] we list 
different predictions for masses of exotic baryons and discuss the N c count- 
ing for the mass splittings. In Sect. E| we repeat the same analysis for the 
decay widths. Summary is given in Sect. |SJ 

2. Soliton models 

Soliton models are often regarded as orthogonal to the quark picture. 
Very often they are generally referred to as Skyrme type models where 
only mesonic degrees are present. In this Section we will demonstrate that 
they are deeply rooted in QCD, take into account quark degrees of freedom 
maybe even in a more complete way than the quark models themselves, and 
that they are fully operative providing predictions of static baryon proper- 
ties, structure functions, skewed and off-forward amplitudes and light-cone 
distribution amplitudes for baryons (for review see e.g. Refs. [T31 [TBI H7\ ). 
not to mention properties of pseudoscalar mesons 18 . That of course does 
not mean that they capture all physics, since — for example — they do 
not posses confinement. They rely on large N c limit and chiral symmetry 
breaking. We shall also make distinction between quark-soliton and Skyrme 
model. 

Let us take as a starting point the chiral Lagrangian density of the form 

C = 4)(i$ - MU Tj [(p\)il) (1) 

which looks like a Dirac Lagrangian density for a massive fermion ip if not 
for matrix U. In fact ip is a 3- vector in flavor space and also in color. Matrix 

JJ15 = e (i/F v )X-0-/ 5 

parameterized by a set of pseudoscalar fields (p has been introduced to re- 
store chiral symmetry given by a global multiplication of the fermion field 
by a phase factor 

i>^e il ' d ^i!. (2) 

Indeed, the term Mtpif) is not invariant under @, however M ipU lb [ipjip is, 
provided we also transform meson fields 

UT*[<p] _> e-^^U^itp] e" ild ^. (3) 
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Since matrix U "lives" in the flavor space, the color indices are here con- 
tracted producing simply an overall factor N c in front of 

Lagrangian Q does not contain kinetic term for meson fields, so <^'s are 
expressed in terms of fermion fields themselves. The kinetic term appears 
only when we integrate out the quark fields. Then the resulting effective ac- 
tion contains only meson fields and can be organized in terms of a derivative 
expansion 



where constants F v and e can be calculated from JQ) with an appropriate 
cut-off. Avz is the Witten Wess-Zumino term which takes into account 
axial anomaly and does not require regularization. Perhaps the most im- 
portant part are the ellipses which encode an infinite set of terms that are 
effectively summed up by the fermionic model of Eq. (^Q). The truncated 
series of Eq. (J3J) is the basis of the Skyrme model. Hence the Skyrme model 
is (a somewhat arbitrary, because it does not include another possible 4 
derivative term) approximation to (Q). 

At this point both models, chiral quark model of Eq. (^Q) and Skyrme 
model of Eq. (jlj) (without the "dots" ) , look like mesonic theories devised to 
describe meson-meson scattering, for example [TJ||. Baryons are introduced 
in two steps, following large N c strategy described by Witten in Refs. |2U| . 
First, one constructs a soliton solution, i.e. solution to the classical equa- 
tions of motion that corresponds to the extended configuration of the meson 
fields, i.e. to matrix U which cannot be expanded in a power series around 
unity. Second, since the classical soliton has no quantum numbers (except 
baryon number, see below), one has to quantize the system. Perhaps this 
quantization procedure, which reduces both models to the nonrelativistic 
quantum system analogous to the symmetric top |211 [1] with two moments 
of inertia I\ 2, makes chiral-soliton models look odd and counterintuitive. 

It is not our purpose to give the full review of the soliton models which 
can be found elsewhere ^lEUE!) especially their connection with QCD, 
to some extent obvious from Eq. (Q), was extensively reviewed in Ref. |22| . 
Here we want to discuss the interconnection of the chiral soliton models and 
quark models. In fact Lagrangian (^Q) is an interacting quark model, despite 
the fact that there are no gluons there. The interaction proceeds through 
the chirally invariant coupling ipU^ip and information about gluons (which 
are integrated out) is encoded in the coupling strength M (constituent quark 
mass) . 




Tr d„U tft, d v Utf + r wz + . - - , (4) 
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For the purpose of illustration it is convenient to use the variational 
approach for the soliton solution [23]. To this end one uses a hedgehog 
ansatz for the static Uq field: 



e in-fP(r) o 
1 



(5) 



where 2x2 matrix in the upper left corner depends on one variational 
function P(r) = P(r/rg) characterized by an effective size rg and being a 
subject to the boundary conditions P(0) = tt and P(oo) = 0. For ro = 
matrix Uq = 1 and the spectrum of the Dirac operator corresponding to (^Q) 
looks like a spectrum of a free fermion of mass M (see the right panel of 
Fig. Once we increase the size ro, the levels rearrange and one distinct 
level "sinks" rapidly into a mass gap (see the left panel of Fig. U) . 
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Fig. 1. Spectrum of the Dirac operator in the presence of the valence level and 
without. The soliton energy is calculated as the regularized difference of two con- 
tributions 

The energy of this interacting fermionic system is given as a sum of the 
valence level and the sea levels (filled levels of the Dirac "sea") calculated 
with respect to the vacuum (ro = 0) configuration as depicted in Fig. ^ 
Stable minimum is achieved for some intermediate soliton size ro = r so \, 
usually of the order of a fermi, as an interplay between decreasing energy 
of the valence level and increasing energy of the Dirac sea. In this case the 
baryon number of the soliton is given simply as the baryon number of the 
valence level. 

An interesting limit can be considered by artificially tuning the size 
of the soliton to ro = |24| . In this limit the valence level goes back 
to the upper edge of the mass gap and the contribution of the sea levels 
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cancels out. Hence the soliton in this limit looks like the constituent quark 
model. Indeed, it has been shown that for ro — > many static observables 
calculated in the soliton model are in agreement with the naive quark model 
predictions. These include qa = 5/3, AS = 1 and /u n /// p = —2/3. 

In the Skyrme model the soliton is constructed purely from the mesonic 
field ©• The baryon number is given as a charge of the conserved topo- 
logical current. Stabilization is achieved by an interplay of the increasing 
energy of the kinetic term and the decreasing energy of the Skyrme term (0J) . 
This reflects the main difference between quark-soliton and Skyrme-soliton. 
Had we included all terms denoted by ellipses in Eq. (^Q) there would have 
been no minimum of energy as a function of ro. 

The quantization on the other hand proceeds in both models almost 
identically [3]. The symmetric top Hamiltonian is supplemented by a con- 
straint coming from the Avz term which selects SU(3) representations that 
contain states of Y = 1. Octet, decuplet, exotic antidecuplet and eikosihep- 
taplet (i.e. 27-plet) are the lowest possible representations satisfying 
this constraint. The difference between the two models is buried in the an- 
alytical form of the expressions for the symmetric top parameters (overall 
mass and moments of inertia, etc.). Some of them are identically zero in the 
Skyrme model, whereas they are nonzero in the quark-soliton model due to 
the valence level contribution. 

Is the tower of representations satisfying constraint Y = 1 infinite? For- 
mally the answer is: yes, but physically: no, since we have to revise as- 
sumptions which led us to the quantization of the soliton. Namely, we have 
assumed rigid rotation which is (classically) very unlikely when the soliton 
angular velocities become large. Two phenomena are expected: deformation 
of the soliton and vibrations. Deformation will lead to instabilities resulting 
in radiation of pions (Goldstone bosons in general). Fast rotating solitons 
will have a cigar-like shape and will lie on linear Regge trajectories |26| . 

As discussed above there is only one representation of given dimension 
in the allowed series of representations selected by the Wess-Zumino con- 
straint. So there is only one (nonexotic) octet, while the quark models in- 
evitably require a cryptoexotic octet together with antidecuplet. Of course, 
the octet is not missing; it has to be of different origin. So far we have 
constructed only rotational states, however, there will be also vibrations. 

Explicit construction of the vibrational states in the Skyrme model (with 
the dilaton field) has been carried on by Weigel in Refs. ;8 . In this approach 
only one mode, namely the "breathing" mode of the soliton was quantized 
and a subsequent mixing with other states was investigated. 
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3. Mass estimates 

In order to estimate the mass of + we have to know two quantities, 
namely the strange moment of inertia, I2 which contributes to antidecuplet- 
octet splitting 

AM W _ 8 = A ( 6 ) 

and the energy shift due to the nonzero strange quark mass. It has been 
observed by Guadagnini [27] that the mass splittings in the Skyrme model 
are well described in terms of 2 parameter effective Hamiltonian 

H' = aD<§+j3Y, 

where a ~ S n N- However (5 = in the minimal Skyrme model and the 
spectrum cannot be well described in the first order perturbation in H' . 
Second order correction, which can be schematically written as [5] 

AE 2 ' h x a 2 (7) 

mimics the nonzero (5. This second order correction cannot be smaller than 
the typical mismatch of the first order AM. Therefore, there is a lower 
bound on AE2, which translates into an upper bound on AMj^_ g . The 
antidecuplet cannot be too heavy. On the other hand AE% cannot be too 
large for consistency reasons, hence too light antidecuplet is also excluded. 
This is how the original prediction Mq+ ~ 1535 MeV was obtained [Hj. The 
updated results of this analysis are given in Table I. 

In the quark-soliton models one chooses different path. Instead of going 
to the second order in perturbative expansion in m s one calculates non- 
leading terms in 1/N C [2SJ. This generates (3 7^ from the beginning and 
the lower bound on I2 does not exist. One can try either to constrain I2 
by identifying some known nucleon resonance with , as it was done in 
Ref. [7], or resort to explicit model calculations which, however, cover rather 
broad range of allowed values 

In the original paper of Diakonov, Petrov and Polyakov [7j the value 
of I2 was fixed by the identification of with A r *(1710) and the equal 
spacing in antidecuplet by adopting the value of 45 MeV for E^n- 

Today we would take a different approach. We would use + rather 
than iV^y to fix the average 10 mass. In a recent paper [2] it has been shown 
that the set of parameters of the symmetry breaking Hamiltonian 

H' = aD^+(3Y + ^=D^Si (8) 

(8) 

(where Dg g are SU(3) Wigner matrices, Y is hypercharge and Si is the 
collective spin operator [2B]) which reproduces well the nonexotic spectra, 
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TABLE I 

Masses of baryons obtained by minimizing the square deviation with respect to 
M so \, Ii I2 for fixed a. 
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Fig. 2. Antidecuplet of SU(3) flavor including <9 + . In Ref. [7] I2 was fixed by 
JV*(1710) and the splitting 8 by fixing S nN . 

as well as the measured mass of the <9 + (1540), can be parametrized as 
follows 1 : 

= 336.4- 12.9^, (5 = -336.4 + 4.3^, 7 = —475.94 + 8.6 . 

(9) 



1 We use here m 3 /(m u + m d ) = 12.9 
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Moreover, the inertia parameters which describe the representation splittings 

Z\M 10 - 8 = J- , AM W _ S = J- (10) 

take the following values (in MeV) 

-^ = 152.4, — = 608.7 - 2.9 17,^ • (11) 

If, furthermore, one imposes additional constraint that = 1860 MeV, 

then E n N = 73 MeV |5| (see also [30] ) in agreement with recent experimental 
estimates [3*T] . 

Hamiltonian (JHJ) introduces mixing between different representations 

\B 8 ) = |8 1/2) 5)+4|T0 1/2 ,B) + 4|27 1/2)J B) , 
|Sio) = |l03/ 2 ,B>+af 7 |27 3/2 , J B> + af 5 |353/ 2 ,5> , 
\Bjo) = |T0 1/2 , B) + df |8 1/25 B) + df 7 |27i /2 ,S> + d§ |35 1/2 , 5) , (12) 

where denotes the state which reduces to the SU(3) representation 1Z 
in the formal limit m s — > 0. The m s dependent (through the linear m s 
dependence of a, and 7) coefficients c^, c?£ and in Eq. (|T2*|) can be 
found e.g. in Ref. [SJ. 

Although the model seems to describe the spectrum of antidecuplet 
rather well (assuming that <9 + and Sjq exist and have masses as discussed 
in the Introduction), we encounter here the first potential problem. Namely, 
the exotic-nonexotic mass splitting © reads in fact 

^2 

whereas 

AM 10 -s = ^- ~ ©(A^ 1 ) • (14) 

This iV c counting is in fact born by experiment 

ZiMio-s - 230 MeV , 4A%_ g ~ 600 MeV , (15) 

however, it poses a serious problem to the validity of the quantization pro- 
cedure for exotic states in large N c limit. Indeed, in the imaginary world of 
very large N c all nonexotic states are degenerate, whereas exotic ones are 
split by a quantity of the order 0(1), similarly to the vibrations with which 
they will mix. This fact although known already in the late 80's have been 
recently revised critically in the literature |121 113j 



epiphany printed on February 2, 2008 



9 



4. Decay widths 

In the decay of + —* NK the kaon momentum in the rest frame of the 
decaying particle 

p K = 267 MeV (16) 
is almost identical to the pion momentum in A decay 

p n = 225 MeV. (17) 

One would, therefore, naively expect that the decay width of (9 + should be 
at least as large as the one of A or even larger, since no suppression com- 
ing from the overlap of the wave functions is expected. Indeed, pentaquark 
states have — naively — so called fall- apart modes. One of the chief pre- 
dictions of the quark-soliton models is that the <9 + width, contrary to the 
naive expectations, is very small [7j. This prediction stimulated experimen- 
tal searches. 

Whilst the mass spectra discussed in the previous section are given as 
systematic expansions in both N c and m s in a theoretically controllable way, 
reliable predictions for the decay widths cannot be organized in a similar 
manner. In fact the decay width is calculated by means of the formula for 
the decay width for B — > B' + ip: 

up to linear order in m s . The "bar" over the amplitude squared denotes av- 
eraging over initial and summing over final spin (and, if explicitly indicated, 
over isospin) . Anticipating linear momentum dependence of the decay am- 
plitude M. we have introduced reduced amplitude A which does not depend 
on the kinematics, i.e. on the meson momentum p^. For the discussion of 
the validity of (|TH|) see 

Soliton models can be used to calculate the matrix element Ai. In order 
to match former normalization j7j we shall define the decay amplitude as 



M B ^B'^ = (B , \6^\B) 

= 3 (B'\ G D vl - G x d ibc D { ^S C - ^D®§i \B) x p% , (19) 



where the sum over the repeated indices is assumed: i = 1,2,3 and 6, c = 
4, ... 7. It is assumed that coupling constants G^i^ can be related to the el- 
ements of the axial current by means of the generalized Goldberger-Treiman 
relations [7j- Explicitly 



2 



_ 3G n n 3 
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For antidecuplet decays (1Z = 10): 



^To _ - Gi - -G 2 , C 1 q+^ n+k — -, (20) 



whereas for decuplet (1Z = 10): 



Gio - G + -G 2 , C^ n+7T - - . (21) 



Hence the suppression of antidecuplet decay width may come only from the 
cancellation between Go,i,2 entering Gjq. Indeed, in the nonrelativistic small 
soliton limit discussed above one can show that G\/Gq = 4/5, G 2 /Gq = 2/5 
and Gjq = 0! Although nonintuitive this cancellation explains the small 
width of antidecuplet as compared to the one of 10 for example. 
One problem concerning this cancellation is that formally 

G ~ 0(iV c 3/2 ) + 0{N]/ 2 ), Gi, 2 ~ 0{N]/ 2 ) (22) 

and it looks as if the cancellation were accidental as it occurs between 
terms of different order in N c . That this is not the case was shown in 
Ref. 33 . Indeed for arbitrary N c antidecuplet 10 = (0,3) generalizes to 
"TO" = (0, ^fcfcS), decuplet "10" = (3, ^) and octet "8" = (1, &f±) [33j, 
and the pertinent Clebsch-Gordan coefficients in fact depend on N c : 



C"To" ~ G o "J - Ci - -jG 2 ■ (23) 



So the subleading Gi-term is enhanced by additional factor of N c and the 
cancellation is consistent with N c counting. Moreover 

m , 3(JV C + 1) / 1 



e+ ^ N + K (Ar c + 3)(AT c + 7) ~ \N C 
,« 10 » _ (JV c -l)(JV c + 5) 
2(N c + l)(N c + 7) 



C A ^ N+7T = n/Ar , lWAr , ^ ~ (24) 



and it looks like the antidecuplet width were suppressed with respect to 
decuplet. Unfortunately the phase space factor spoils this counting. 
Indeed, because of (|T3*|) and (|Tl*|) 

Pn ~ O (^-] , pk ~ 0(1) (25) 



and consequently 



ij-.v+^-^jij]. r„:_. v+K -~(o ( i) (20) 
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Fig. 3. The correction factors i?( mlx ) due to SU(3)-breaking representation mixing 
for the decays discussed in the text, as functions of the parameter p = Gi/Gq. 

in the chiral limit. This N c counting contradicts experimental findings which 
suggest r e +^ N+K < r A ->N+n- 

A few comments are in order. Firstly, let us note that in Nature neither 
7T nor K mesons are massless and both p n and px are of the same order 
of 230 MeV (|Tfi)l . ([T7|) and the scaling does not hold. So for m n ^ 
and rriK ^ both meson momenta scale as however A does not decay, 
because in the large N c limit it is degenerate with nucleon and cannot emit 
a massive particle, whereas + does decay. It is an instructive example how 
subtle is an interplay of theoretical limits N c — > oo and m q — > 0. Secondly, 
((^51) holds only if the cancellation Gjq = is not exact. Let us suppose 

1/2 

that the leading N c power cancels in such a way that G w ~ 0{Nc' ) rather 

than C(iV c 3/2 ). That would make r +^ N+K ~ C(1/7V C 2 ), i.e. of the same 
order as r^N +7T . Finally let us remark that there is further suppression 
of Fq+_> n+k coming from the mixing (|12fl . This is illustrated in Fig. 01 
where we plot multiplicative correction factor i?( mix ) as a function of the 
parameter p = G\/Gq. Phenomenological value of p ~ 0.5 (Hj. 

5. Summary 

The solitonic approach to baryons is very successful, as it describes not 
only spectra and other static properties, but also structure functions, skewed 
and nonforward parton distributions and also ligth-cone distribution ampli- 
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tudes. However, it relies on many approximations. Firstly, it is based on an 
ansatz and as such one must check the self-consistency of the approach. Ev- 
idently, following arguments by Witten [20], solitonic solutions are justified 
only for large N c . Indeed, looking e.g. at Eq. ((J) a calculation of bary- 
onic properties requires functional integration over both bosonic ((/?) and 
fermionic (ip) degrees of freedom. In order to apply the stationary phase 
approximation, as it is commonly done in the soliton models, one has to 
omit the bosonic functional integral, using instead the background bosonic 
field that minimizes the effective action of the system. This is only justified 
for large iV c , where the bosonic fluctuations are suppressed. Secondly, soli- 
ton quantization proceeds by quantizing the rotations in space and flavor 
space. To this end one assumes the rotational motion to be adiabatic. This 
means that angular velocities go like J /I ~ 1/N C yielding frequencies (and 
hence excitation energies) of order 0(N~ 1 ). This, in turn, implies a Born- 
Oppenheimer separation of the slow collective rotational motion from the 
faster modes associated with vibrations. Because of this scale separation 
the collective rotational modes can be quantized separately from the intrin- 
sic vibrations. While this procedure has been applied with great success to 
many properties of the nonexotic baryons it has been criticized as far as 
exotic multiplets are concerned. 

The question j!21 IP3] here is whether the rigid-rotor type semiclassical 
projection can be applied to exotic states. The fact that the standard 
semiclassical quantization gave excitation energies of the order 0(N®) for 
exotic states (fT3)) means that the approach is not justified for such states 
unless further arguments can be invoked. In contrast, in view of Eq. 1)141) . 
rigid-rotor quantization is certainly justified for the non-exotic states. 

Diakonov and Petrov .35. have argued that while it is true that at large 
N c rotational excitations are comparable to vibrational or radial excitations 
of baryons, both non-exotic and exotic, the corrections due to the coupling 
between rotations and vibrations die out as 1/N C . The collective rotational 
quantization description fails only when the exoticness, i.e. the number of 
qq valence pairs needed to construct the quantum numbers of a given state, 
becomes comparable to N c . The newly discovered + baryon belongs to the 
antidecuplet of exoticness = 1. The larger N c , the more accurate would be 
its description as a rotational state of a chiral soliton. Diakonov and Petrov 
|35| support their estimates by considering a simple model consisting of a 
charged particle in the field of a monopole. However, if one generalizes 
the model by considering two charged particles interacting by a harmonic 
potential and moving in the field of a monopole, the coupling of rotational 
and vibrational degrees of freedom of this system is by no means vanishing 
but strong [T5] . 

As we have discussed in Sect. 0] large- N c arguments apply also to the 
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width of the baryons considered, because if the approach is justified, then 
at a formal level the width must approach zero at large N c . Of course, for 
non-exotic states such as the decuplet, this is true. The reason is simply 
phase space (|25|). Unfortunately, as shown recently in Ref. [331 the width 
of the + as calculated via the standard collective rotational approach is 
of the order N® in the chiral limit. This demonstrates that the procedure 
is not self consistent on the basis of pure large- iV c arguments. Thus, if the 
width of the + is really small, as the experiments indicate, there must be 
particular dynamical reasons for the smallness, which exist on top of what 
is required for the validity of the large- N c expansion alone. In this context 
the cancellation leading to Gjq = in the small soliton limit is of particular 
importance. 

For completeness one should also mention another approach to the quan- 
tization of chiral solitons based on the assumption that SU(3) symmetry is 
strongly broken [Tl|. This approach, known as a bound-state approach, was 
recently applied to + by Klebanov et al. These authors reconsider 

the relationship between the SU(3) rigid-rotator and the bound-state ap- 
proach to strangeness in the chiral soliton models. For non-exotic S = — 1 
baryons the bound-state approach matches for small kaon mass nix onto 
the rigid-rotator approach, and the bound-state mode turns into the rota- 
tor zero-mode. However, for small mx, there are no S = +1 kaon bound 
states or resonances in the spectrum (unless txik = 0). This shows that 
for large N c and small (but non-zero) itik the exotic state is an artifact 
of the rigid-rotator approach. An S = +1 near-threshold state with the 
quantum numbers of the + pentaquark comes into existence only when 
SU(3) breaking is sufficiently strong or vector mesons are introduced |3b| . 
Therefore, one argues that pentaquarks are not generic predictions of the 
chiral soliton models. 

The present work is supported by the Polish State Committee for Sci- 
entific Research (KBN) under grant 2 P03B 043 24. 
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